Spin-Peierls transition in TiP04 
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We investigated the magnetic and structural properties of the quasi-one dimensional -quantum 
chain system TiP04 (J 965 K) by magnetic susceptibility, heat capacity, ESR, x-ray diffraction, 
NMR measurements, and by density functional calculations. TiP04 undergoes two magnetostruc- 
tural phase transitions, one at 111 K and the other at 74 K. Below 74 K, NMR detects two different 
^^P signals and the magnetic susceptibility vanishes, while DFT calculations evidence a bond alter- 
nation of the Ti. . . Ti distances within each chain. Thus, the 74 K phase transition is a spin-Peierls 
transition which evolves from an incommensurate phase existing between 111 K and 74 K. 
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The discovery of high- Tc superconductivity in two- 
dimensional oxocuprates has stimulated broad interest 
in the properties of low-dimensional quantum 3=1/ 2 an- 
tiferromagnets. The complex interplay between spin, 
charge, orbital and lattice degrees of freedom in low- 
dimensional systems with pronounced quantum fluctu- 
ation renders a plethora of complex and unusual ground 
states. [H, [3 

Most of the prominent examples of low-dimensional 
quantum antiferromagnets with exotic ground states con- 
tain Cu^+ {3dP, S = 1/2) ions with one hole present 
in the Cg orbitals.[3|-[5| Compounds of early transition- 
metal elements with one electron in the d shell are less 
frequently investigated. With a 3d^ electron in a high- 
symmetry or slightly distorted octahedral environment, 
the orbital degeneracy of the t2g states opens new degrees 
of freedom, with the possibility of low-energy orbital ex- 
citations and the interesting scenario of destabilization 
of coherent spin/orbital ordering by quantum fluctua- 
tions. 3d^ systems can be easily realized in compounds 
containing, e.g., Ti^+ or cations. A paramount ex- 
ample is the vanadium ladder compound a'-NaV2 05.[6| 
At high temperatures a'-NaV205 contains mixed- valent 
vanadium cations with one electron occupying an or- 
bital confined to the rungs of the ladder, hence constitut- 
ing a quarter-filled ladder system. [7] Below ~ 34 K a'- 
NaV205 undergoes charge ordering [8] leading to a spin- 
gap of ~ 100 K indicated by a rapid drop of the mag- 
netic susceptibility. That was initially ascribed to a spin- 
Peierls transition. [gI, 111 Other prominent low-dimensional 
3d^ systems that have recently attracted much atten- 
tion are the Mott insulators TiOX (X=C1, Br).[io|-[ii| 
These compounds crystallizing in the FeOCl-type struc- 
ture, consisting of Ti - O - X layers made up of Ti04Cl2 
octahedra. These layers are stacked with van der Waals 
interactions between them. [15, 16] The magnetic suscep- 
tibility of TiOCl reveals several unusual features, which 
led to the proposal that TiOCl may be a manifestation of 



a resonating valence-bond solid. [10| Subsequently, Seidel 
et al. demonstrated that the high temperature suscepti- 
bility fits very well to a 8=1/ 2 Heisenberg chain model 
with nearest-neighbor (nn) antiferromagnetic (afm) spin- 
exchange (SE) interaction of ~660 K. In view of these 
findings and their LDA + U electronic structure calcu- 
lations, Seidel et al. concluded that TiOCl is an exam- 
ple of a Heisenberg chain that undergoes a spin-Peierls 
transition at 67 K.[ll| Subsequent low temperature x-ray 
structure determination showed a slight dimerization of 
the Ti. . . Ti distances along the h direction, [l^ A further 
anomaly was detected at Tc2 = 95 K, which signals a 
transition into an incommensurate phase with a slight 
monoclinic distortion of the lattice. It was initially be- 
lieved and it appears to be now generally accepted to be 
of continuous or higher order, however, there was a later 
claim that it was first order. [13, [H-Ell 

Here, we report the magnetic and structural proper- 
ties of TiP04, which contains Ti^+ cations and displays 
two magneto-structural phase transitions reminiscent of 
those in TiOX. In contrast to TiOCl, however, TiP04 is 
a structurally one-dimensional compound crystallizing in 
the CrV04 structure-type (SG: Cmcm) (see inset in Fig. 
[T]).[22] The Ti+^ ions, carrying 8=1/ 2 entities, are sub- 
ject to axially compressed TiOe octahedra. These share 
their edges to form corrugated Ti04 ribbon chains along 
the c-axis with a buckling angle of 156.927(4)° in the a-c 
plane. The Ti^+...Ti^+ distance at room temperature 
(RT) amounts to 3.1745(10)A with a Ti3+-02--Ti3+ Z 
of 95.484(5)°. [m The Ti04 ribbon chains are intercon- 
nected by sharing corners with distorted PO4 tetrahedra. 

At high temperature the magnetic susceptibility of a 
polycrystalline sample (Fig. (TJ main panel) is character- 
ized by a broad maximum centered at ~ 625 K, indicating 
short range afm correlations. After correction for a tem- 
perature independent offset to the susceptibility arising 
from diamagnetic contributions of the closed shells and 
van Vleck terms, the high temperature magnetic suscep- 
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FIG. 1: (Color online) (o) Molar magnetic susceptibility, Xm, 
of TiP04 measured in a field of IT. The (red) line is a fit 
to a Heisenberg chain with uniform nn afm SE interaction, 
see text. Insets: (a) Crystal structure of TiP04, where green 
and blue polyhedral represent the TiOe octahedra and PO4 
tetrahedra, respectively, (b, c) Xm in the region of the phase 
transitions, (d, e) Heat capacity, Cp, in the region of the 
anomalies, where the red circles and blue triangles refer to 
the heating and cooling data, respectively. 
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FIG. 2: (Color online) ^^P MAS NMR spectra for TiP04 
(temperatures indicated). The (orange) hashed area high- 
lights the incommensurate continuum. Insets: (a) Peak po- 
sitions and/or boundary edges versus temperature, (b) The 
41.1 K spectrum. 



tibility can be described very well by a 5=1/2 Heisen- 
berg chain with a uniform nn afm SE interaction [24] of 
965(10) K and a ^/-factor of 1.94(3). 

Below ~ 120 K the susceptibility reveals two subse- 
quent magnetic phase transitions, indicated by two rapid 
drops of the susceptibility at 111(1) K and 74(0.5) K. Fi- 
nally, at lowest temperatures the susceptibility levels off 
to a value of 75(10) xlO~^ cm^/mol. At very low tem- 



perature a slight increase is seen, which we ascribe to a 
Curie tail due to '-«^70 ppm of a free 5=1/2 spin entities. 
The anomaly at 74 K shows a thermal hysteresis with 
a temperature difference of ~ 50 mK between the heat- 
ing and cooling traces while heating/cooling cycles gave 
identical susceptibilities for the 111 K anomaly (see Fig. 
[1] inset (b, c)). 

Heat capacities collected on crystals exhibit two A- 
type anomalies at 110.9(0.6) K and 74.1(0.3) K, with the 
lower temperature anomaly also showing a thermal hys- 
teresis, while again no hysteresis is seen for the higher 
temperature anomaly (see Fig. [1] inset (d, e)). Angu- 
lar and temperature dependent Electron Spin Resonance 
measurements (ESR) on single crystals revealed a sin- 
gle Lorentzian resonance line (^-factor 1.93 - 1.95) and a 
linewidth decreasing linearly with temperature (50 Oe < 
FWHM < 300 Oe) consistent with earlier findings. [H,!!!! 
The integrated intensity of the ESR line mimics the tem- 
perature dependence of the magnetic susceptibility and 
drops to zero below 74 K. 

The crystal structure of a very high quality single 
crystal of TiP04 was determined by x-ray single crys- 
tal diffraction measurements at various temperatures be- 
tween 293 K and 90 K. Down to 90 K the structure was 
found to be identical to that reported by Glaum et al.^ 
except for small changes of the lattice parameters and the 
general atomic positions. [23| As the temperature is low- 
ered from RT to ~120 K, the lattice parameters a and 
b increase slightly but the c parameter decreases, such 
that the cell volume remains almost constant. The resid- 
ual electron density, i.e. the measured electron density 
minus the calculated electron density (from the super- 
position of spherical atom densities) shows considerable 
residuals located within the a - c plane next to the Ti 
atoms, at an interstitial position in the ribbon chain, and 
at a position bisecting the O - Ti - O angle perpendicular 
to the ribbon chains. Upon cooling there is a gradual mi- 
gration of the density away from the bisecting position 
into the interstitial position within the ribbon chains. 
This change is complemented by a reduction of the dis- 
tance in the Ti - O - Ti, nn super-exchange pathway, and 
an increase of the intrachain O - Ti - O angle. Evidence 
for a structural change was not found in this tempera- 
ture range, possibly due to the dynamic character of the 
intermediate phase as observed by NMR (see below). 

Magic angle spinning (MAS) ^^P nuclear magnetic res- 
onance (NMR) spectra (center field ^ 8.5 T) operating 
at a spinning frequency of ~ 25 kHz, were collected on 
a polycrystalline sample between ~ 35 K and RT. The 
spectra are displayed versus temperature in the main 
panel of Fig. [2j Above ~ 140 K we observe a sin- 
gle ^^P symmetric NMR line accompanied by two sets 
of very weak symmetrically placed spinning sidebands. 
Near 116 K the line becomes asymmetric and below 111 K 
it broadens into an asymmetric continuum limited by 
two boundary peaks. With decreasing temperature the 
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continuum expands and its intensity decreases. Towards 
~ 76 K the continuum finally washes out, whereupon its 
lower boundary grows into two symmetric lines indicating 
the occurrence of two different P atom environments (see 
inset Fig. [2](b)). The peak positions and/or boundary 
edges are shown versus temperature in Fig. [2] (a). There 
are similarities between NMR measurements of TiP04, 
reported here, and TiOX, reported by Saha et a/.. [26] 

We now probe the SE interactions of TiP04 by per- 
forming mapping analysis based on density functional 
calculations. |27| We consider the nn and next-nearest 
neighbor (nnn) intrachain SE interactions Ji and J2, re- 
spectively, as well as the interchain SE interaction J3 (see 
Fig. [3j). To evaluate Ji - J3, we determine the relative en- 
ergies of the four ordered spin states, FM, AFl, AF2, and 
AF3 shown in Fig. [3l by density functional theory (DFT) 
electronic band structure calculations. Our DFT calcu- 
lations employed the Vienna ab initio simulation pack- 
age (Ref. I28I [29I ) with the projected augmented- wave 
method, the generalized gradient approximation (GGA) 
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for the exchange and the correlation functional. [30| We 
used a plane-wave cut-off energy of 400 eV, a set of 56 
k-point irreducible Brillouin zone, and the threshold of 
10 ~^ eV for the self-consistent-field convergence of the 
total electronic energy. To account for the electron corre- 
lation associated to the Ti 3d state, we performed GGA 
plus onsite repulsion (GGA+/7) calculations (Ref. [3l|) 
with an effective Ueff = U - J = 2 eV and 3 eV on Ti. 
The relative energies, per four formula units (FUs), of 
the four ordered spin states are summarized in Fig. [3l 

The total SE energies of the four ordered spin states 
can be expressed in terms of a Heisenberg spin Hamil- 



tonian, H 



y~] J ij Si Sj , 



where Ji\ is the SE interaction 



between the spins Si and Sj on the spin sites i and j, re- 
spectively. By applying the energy expressions obtained 
for spin dimers with N unpaired spins per spin site (in 
the present case, N = 1),[32], the total SE energies for 
the four configurations, per four FUs, are given in Fig. 

n 

Thus, by mapping the relative energies of the four or- 
dered spin configurations given in terms of the SE pa- 
rameters (see Fig. [3]) onto the corresponding relative 
energies obtained from the GGA+/7 calculations, we ob- 
tained the values for the SE parameters Ji - J3 (see Table 
[Tj) . [13, [111 The results of our DFT calculations are in very 
good quantitative agreement with our experimental find- 
ings indicating a very large nn intrachain afm SE inter- 
action. The nnn intrachain SE interaction is almost two 
orders of magnitude smaller, the interchain interaction 
J3 amounts to 2% of Ji. 

The low temperature MAS data of TiP04 prove a 
non-magnetic ground state with two distinct P atomic 
environments, as evidenced especially by the low tem- 
perature spectra. The chemical shifts of the ^^P lines 
amount to ~ -20 ppm in good agreement with what 
has been found for other diamagnetic orthophosphates. 



FM (0,0) 
Efm=(-Ji-J2-J3)/4 



AFl (-42.5,-32.5) 

Eafi=(+Ji-J2-J3)/4 



AF2 (-21.3,-16.3) 

EaF2=(+J2-J3)/4 



AF3 (-43.5,-33.0) 

EaF3=(+Ji-J2 + J3)/4 



FIG. 3: (Color online) The four ordered spin configurations, 
FM, AFl, AF2 and AF3, used to extract the values of Ji, 
J2 and J3, where only the Ti^^ ions are shown for simplicity. 
The up- and down-spin Ti^^ sites are indicated by different 
colors. The numbers in parenthesis (from left to right) repre- 
sent the relative energies in meV per four FUs obtained from 
GGA+t/ calculations with Uef£ = 2 and 3 eV, respectively. 
The expression of the total SE energy per four FUs is also 
given. 



J^ 


U = 2eV 


U - 3eV 


Jl 


-988 


-751 


J2 


-1.4 


+0.7 


Js 


-20 


-15 



TABLE I: Values of the SE parameters Ji - J3 derived from 
the mapping analysis (in K). 



proving the non-magnetic character of the ground state 
of TiP04. 34, 35I We ascribe the 74 K phase transition in 
TiP04 to a spin-Peierls transition with the Ti. . . Ti bond 
alternation within the Ti chains. 

To probe the low temperature crystal structure of 
TiP04, we considered the subgroups Amm2 and Pmmn 
of the RT space group Cmcm. By GGA calculations, 
we optimized the structures of TiP04 starting with the 
initial settings described by Cmcm, Amm2 and Pmmn 
without symmetry constraints in order to allow the atom 
positions to relax freely (with a set of 28 k-point irre- 
ducible Brillouin zone, and the thresholds of 10~^ eV 
and 0.001 eV/A for the self-consistent-field convergence 
of the total electronic energy and force, respectively). 
The lowest-energy structure found was obtained start- 
ing from the Pmmn initial setting, and was lower in 
energy by ~ 48 meV per formula unit than was the 
Cmcm structure, and by ~ 32 meV per formula unit 
than the Amm2 initial structure. The structure relaxed 
from the Pmmn initial setting shows a dimerization in 
the Ti chains with alternating Ti. . . Ti distances of ~2.9 
and ~3.5 A, which is comparable in magnitude to that 
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observed in TiOCl.[l3 This structure also has two dif- 
ferent environments for the P atoms within the PO4 
units, which is consistent with our MAS NMR spec- 
tra. From model calculations using the relaxed crystal 
structure data, we expect very weak superstructure re- 
flections, which were not resolved in an early neutron 
powder diffraction experiment. (36| 

The incommensurate phase seen between ~ 111 K and 
~ 74 K is similar to that found for TiOX, where it has 
been attributed to a frustration between the spin-Peierls 
pairing and an elastic interchain coupling. Analogous 
arguments may apply for TiP04 as well. In view of the 
small interchain coupling the incommensurate phase be- 
tween 111 K and 74 K could also be ascribed to a dynamic 
equilibrium between short-range-ordered dimerized seg- 
ments. 

In conclusion, our magnetic susceptibility, heat capac- 
ity, ESR and ^^P MAS NMR measurements supported 
by our density functional calculations show that TiP04 
undergoes a spin-Peierls transition at 74.1(0.3) K, which 
is preceded by an incommensurate phase extending up to 
^ 111 K. The thermal hysteretic behavior of these transi- 
tions is consistent with the pattern of discontinuous and 
continuous transitions seen for TiOCl and TiOBr. At 
high temperatures the magnetic susceptibility of TiP04 
is described by a S=l/2 Heisenberg antiferromagnetic 
chain with an unprecedented nn SE of ~ 1000 K. 
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